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We study forecasts for the accuracy of the determination of cosmological parameters from future 
large scale photometric surveys obtained using the full shape of the 2-point galaxy angular correlation 
function. The effects of linear redshift-space distortion, photometric redshift gaussian errors, galaxy 
bias and non-linearities in the power spectrum are included on our analysis. The Fisher information 
matrix is constructed with the full covariance matrix, including the correlation between nearby 
redshift shells arising from the photometric redshift error. We show that under some reasonable 
assumptions, a survey such as the imminent Dark Energy Survey should be able to constrain the 
dark energy equation of state parameter w and the cold dark matter density ficdm with a precison 
of the order of 20% and 13% respectively from the full shape of the angular correlation function 
alone. When combined with priors from other observations the precision in the determination of 
these parameters improve to 8% and 4% respectively. 

PACS numbers: 95.36.+X, 98.65.Dx, 98.80.Es 

I. INTRODUCTION 

The clustering of large-scale structure of the universe is an invaluable source of information about the fundamental 
parameters that determine a given cosmological model. In the standard cosmological model, its origin is related 
to quantum fluctuations of the inflaton field during the inflationary period, and its evolution is determined by the 
different components of the universe that contribute to the total energy-momentum tensor, such as cold dark matter, 
dark energy, baryons and relativistic particles. 

Until recently, cosmological constraints from galaxy distribution were obtained with spectroscopic catalogues mostly 
from the 2dF m and SDSS projects. These projects raised the study of the galaxy clustering in 3D space to a new 
level, with a main focus on the 2-point correlation function, both in configuration and Fourier space. In particular, 
they resulted in the first observation of a feature in the correlation of the matter ditribution, arising from the so-called 
baryon acoustic oscillations (BAO) that occurred at an early stage of the evolution of the universe 

In the near future, projects such as the WiggleZ Dark Energy Survey ^ and the Baryon Oscillation Spectroscopic 
Survey^ (BOSS) will continue improving the spectroscopic catalogues. Future projects, both ground based such as 
BigBOSS'^ and space based such as Euclid'* will also rely on spectroscopic measurements of objects. 

In order to collect a larger number of galaxies more efficiently, there are a number of planned galaxy surveys that 
will estimate the galaxy redshift from broad-band photometry measurements (photometric redshift or photo-z) with 
the aim to measure more than hundreds of millions up to a few billions of galaxies in a larger volume compared to 
spectroscopic surveys. The idea is to have a trade-off between precise spectroscopic redshifts for a relatively small 
number of galaxies and less precise photometric redshifts for a larger number of objects. For instance, the MegaZ- 
LRG photometric catalogue based on the SDSS-II DR7 has become recently avaliable H, [H]. The measurement of the 
angular correlation function (ACF) for a 1.5 million Luminous Red Galaxies (LRG) with photometric redshifts was 
perfomed in whereas their angular power spectra (APS) was measured in jsj. 

The Dark Energy Survey^ (DES) is a large dedicated photometric survey that will soon start observations. It is 



^ See http://wigglez.swin.edu.au 

2 See http://cosmology.lbl.gov/BOSS 

3 See http://bigboss.lbl.gov 

* See http://sci.esa.int/ouclid 

^ See http://www.darkenergysurvey.org 
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expected to measure 0(300 million) galaxies over an area of 5000 square degrees of the southern sky using 525 nights 
in the Blanco 4-meter telescope in Chile. The main goal of DES is the measurement of the dark energy equation of 
state parameter w up to rcdshift z ^ 1.4. This measurement will be performed using four complementary methods: 
distances from supernovae, large scale structure from the galaxy distribution, number count of galaxy cluster and 
shear from weak leasing measurements 0. 

The purpose of this paper is to determine forecasts for the precision of the determination cosmological parameters 
from future photometric measurements of the large scale structure of the universe, such as what will be obtained by 
DES; using for the first time the full shape of the galaxy correlation function. The use of the angular correlation 
function in a given redshift shell as an observable to infer cosmological parameters can be more appropriate than the 
use of the spatial correlation function in the presence of errors inherent in the photometric measurements of redshifts. 
We use the full shape of the angular correlation function and its covariance, both at different angles and in different 
redshift shells, to perform a Fisher matrix analysis of the sensitivity to cosmological parameters. Hence, we do not 
need to assume any parametric form of the ACF, as done in an interesting recent study 0|. Also our analysis is 
complementary to the study of the angular power spectrum since both techniques have advantages and disadvantages. 
The main disadvantage of the ACF is that the covariance matrix is highly degenerate due to large correlations among 
different angular scales whereas for the power spectrum errors in a given angular scale may propagate into several 
peaks. 

In this paper we model the ACF with 6 cosmological parameters, assuming a spatially flat universe: the dark 
matter density parameter flcdrm the baryon density parameter rif,, the Hubble parameter h, the primordial index of 
scalar perturbations ris , the normalization of perturbations as and the equation of state parameter for dark energy 
w. In addition we allow for different bias between dark matter and galaxies in each of the redshift shells. We include 
in the modelling the effects of linear redshift-space distortion, photometric redshift gaussian errors, galaxy bias and 
non-linearities in the power spectrum. The Fisher information matrix is constructed with the full covariance matrix 
from the model. In particular, we take into account the correlation between nearby redshift shells arising from the 
photometric redshift error. We then discuss various forecasts for the cosmological parameters in different scenarios. 

This paper is organized as follows. In the next Section we present the model for the angular correlation function. 
We include the effects of nonlinearities in the power spectrum, redshift distortions from peculiar velocities of galaxies 
and the photo-z errors. In order to develop some intuition, we discuss in the third Section how the cosmological 
parameters affect the ACF. Section IV is devoted to the modelling of the covariance matrix for the measurement 
of the ACF in different redshift shells with different angular binnings. The Fisher matrix approach is discussed in 
Section V. Our main results arc presented in Section VI and Section VII provides a summary and our conclusions. 



II. MODELLING THE ANGULAR CORRELATION FUNCTION 



The angular correlation function uj{9) is related to the two-point spatial correlation function in redshift space 

by 

iuid)^ dzifizi) dz2f{z2)i^'Hr{zi,Z2,d)). (1) 
Jo Jo 

The function f{z) is determined by the selection function of the survey (/>(2;), the dark matter to luminous bias factor 
b{z) and the linear growth function D{z) (normalized to D{z = 0) = 1) as f{z) = (f){z)h{z)D{z). The radial comoving 
distance r{zi, Z2t9) is the distance between two galaxies at redshifts zi and Z2 separated by an angle 9 and in a flat 
cosmology (which we will assume here) it is given by the relation, 

r = Vx'(^i) + X^{z2) - 2x{zi)x{z2)cos9, (2) 
where x{^i) is the comoving distance of the object i to us (hereafter we use units with c = 1): 

and H{z) is the usual Hubble function determined by the composition of the universe. 

The redshift-space spatial correlation function corrected from the effects of redshift distortions arising from the 
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peculiar velocities of galaxies is given by (in the plane-parallel approximation) 0, [l^ 



l + i[/3(zi) + /3(z2)] + ^/3(^i)/3(^2) 
3 5 



^[/3(zi)+/3(z2)] + ^/?(zi)/?(z2) 



eo(r)Po(M) 
6(0^2(m) 



^/3(^i)/3(^2) 



?4(r)F4(/i). 



(4) 



Here the Pi{fJ.) are the usual Legendre polynomials as a function oi fi ~ d ■ f (cosine of angle between the line of 
sight d and r) and /3(z) = g{z)/b{z) with (7(2) = dlnZ^/dlna. The correlation multipoles are related to the matter 
power-spectrum Pm{k) through: 



6(0 



27r2 



dkk^Pm{k)jiikr) 



The selection function (j){z) is normalized such that: 



dz (j){z) = N 



(5) 



(6) 



where N is the total number of objects per unit solid angle of the survey. One usually slices the survey into n redshift 
bins i in such a way that: 



</.(z) = ^0,(2), 



(7) 



with 



(/)i(z) = n{z)Wi{z), 



(8) 



where n{z) is the number density of galaxies per unit solid angle and per unit redshift and Wi{z) is a window function 
that selects the i-th redshift bin. In general, n{z) is a function of the limiting magnitude of the survey. We will adopt 
a function of the form uM 



n{z) oc (z/z)^e-(^/^")' ' 
with z = 0.5 being the median redshift of the survey and 

M/.(z) = e(z-zne(z'-''-z), 



(9) 



(10) 



where Q{x) is the Heaviside (or step) function. 

However, in surveys where only photometric redshifts are available one should incorporate into the selection function 
the probability P (z^'^jz) of obtaining a true redshift z given that a photometric redshift z^'^ is measured p"3 - [l5j : 



n{z) 



n{z) / dzP" Wi (zP")P(zP''|z) 



dz""' P{z'"^\z). 



(11) 



Of course in the case P{zp^\z) = <5(zph — z), one recovers Eq.®. 

Usually the probability function for a spectroscopically calibrated galaxies is written as a gaussian distribution: 



Pg{z^''\z) 



2'K(7: 



■ exp 



(z - z"" - Zi,.,,)' 

2a? 



(12) 



with the error given by 



a^(z) =5<,(l + z). 



(13) 
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When not mentioned, we set = 0.03 which is expected for LRGs samples [1, [T^. We will also consider here an 
unbiased photo-z and set Zbias = 0. 

Another effect that must be taken into account is the decrease in the power of the ACF due to the nonlinear 



gravitational clustering [17[ • We will follow Crocce et al [13| and model this effect phenomenologically by introducing 



a nonlinear power spectrum PNL{k) and substituting (neglecting the so-called mode-coupling term): 

Pm{k) ^ PjvL(fc) = P™(fc)exp [~r%LeD\z)/2] (14) 

with TjvL ~ 7 Mpc h""'^. Even though this is a simple-minded approach, it is robust depending on the scale someone 
is probing. Crocce et al [l^ showed that this approach is in good agreement with simulations above 40/i~^Mpc, 
therefore our analysis will be valid and applied above this scale. 



III. COSMOLOGICAL INFORMATION IN THE ANGULAR CORRELATION FUNCTION 

Having a model for the angular correlation function one can study the dependencies on the cosmological parameters. 
The model is specified by the usual cosmological parameters. Throughout this paper we assume as fiducial cosmological 
model a flat ACDM universe with parameters as determined by WMAP7 dark matter density parameter ficdm = 
0.222, baryon density parameter fif, = 0.0449, Hubble parameter h = 0.71, primordial index of scalar perturbations 
Us = 0.963, and normalization of perturbations CTg = 0.801. In addition, we set r^vL = 6.6 Mpc h~^, S^^ = 0.03 and 
b = 2^. 

The cosmological parameters can affect the angular correlation function in four different ways: in the primordial 
matter power-spectrum Pm{k), the growth function D{z), the linear redshift-space distortion /3(z) and in the definition 
of comoving distances (geometry). 

The primordial power spectrum is characterized by the parameters ri„i, /i, rig and erg (the dark energy equation 
of state parameter w has a negligible impact on the primordial power spectrum in most cases). The growth function 
depends on r2,„ and w. The redshift space distortion parameter depends mainly on Q,m, and the comoving distances 
are determined by and w. 

It is well known that there are several degeneracies among the parameters and functions, such as bias, ag and the 
growth function. We will show below that the analysis of multiple redshift shells including redshift-space distortion 
can ameliorate this problem [20j . 

In Fig.(IT]) we show how the angular correlation function changes with w and ^cdm around our fiducial cosmology 
and for the redshift shell 1.0 < z < 1.05. From the top panel one can see that the ACF changes mostly its amplitude 
and the location of the BAO feature. As stated above, the primordial power spectrum does not change appreciably 
with w, but the growth factor D{z) and the comoving distance are affected by this parameter. The contribution from 
the growth factor can be easily estimated. Since for this sheh D[w = —1) = 0.618 and D['w = —0.85) = 0.630, the 
growth function can account for an increase of ~ 3.9% in the amplitude. Hence it fall short of explaining the ^ 14% 
increase of the amplitude at intermediate scales. Furthermore, the increase in D{z) can not explain the shift in the 
position of the BAO feature. 

It is noticeable that the BAO feature shifts towards larger angular scale when going from the fiducial cosmology 
w = —1 to w = —0.85, while when w ~ —1.15 it shifts to smaller angular scales. The explanation comes from the fact 
that one changes the geometry when changing w, and distances have different values, changing the "standard ruler". 
When converting the spatial correlation function to find the ACF, the projection depends on the thickeness of the shell 
in comoving distance. If the shell is wider (in comoving distance) the ACF will have less power and the projection 
offset will be higher. For the fiducial cosmology the thickeness of the shell 1.0 < z < 1.05 when converted to comoving 
distance is Ax — 87.4/i~^Mpc, while for w — —0.85 we have Ax ~ 83.6/i~^Mpc, resulting in a modification of ~ 5% 
in the thickness of the comovin g sh ell. This fact also explains why the BAO feature shifts to smaller angles from the 
fiducial value to w = —0.85 (sl. l2lj. It is interesting that such a small change in the shell thickness can generate a 
non-negligible modification in the ACF shape and amplitude. 

The parameter ^cdm has an impact in all four effects listed above and its impact on the ACF is shown in the 
bottom panel of Fig.(IT]). However, most of the cosmological information for ^cdm will come from the primordial 
power spectrum, with modifications in the fraction of baryons /(, — Qb/flm and the product ^mh. 

A complication in the simple effects illustrated above is the presence of photo-z errors, which can mimic some of the 
effects of changing the cosmological parameters, in addition to redshift space distortions, nonlinearities and galaxy 
bias. All these effects are of course taken into account in our results. 
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FIG. 1: Angular correlation function as a function of the angular scale for different values of the dark energy equation of state 
parameter w (top panel) and the cold dark matter density parameter Q.cdm (bottom panel) for the redshift shell 1.00 < z < 1.05 
with photo-z dispersion ~ 0.03. The vertical dotted line indicates the peak position of 6^ x ijj{6) in the fiducial cosmology 
as stated in the text. 



IV. MODELLING THE ERRORS IN THE CALCULATION OF THE ACF 



In order to derive meaningful cosmological constraints or even determine the reach of a given survey, one has to 
understand the errors involved in the measurements of the ACF. More precisely, one must model the full covariance 
matrix for measurements of the ACF at different angles and in different shells. 



A. Covariance matrix for one shell 

Wc will follow closely the works of Blake el al [l^ and Crocce et al [l^ in our modelling. The modelling will 
include effects from photo-z errors, redshift distortion, partial sky coverage and shot noise. For this analysis it is more 
convenient to work with the angular power spectrum (a derivation of the ACF covariance matrix in configuration 
space can be found in Cohn [22|). We define the projected density fluctuations onto the sky in a particular direction 
h, cr(h) as 



oo / 



CT(n) = / dz (j){z) Sin,z) (15) 

and decompose it in spherical harmonics: 

= X! X! O'hnY'hnifl). (16) 
1=0 m=-l 

The angular correlation function can be written as: 

uj{e) ^ {a{n)a(n + ej) (17) 
= ^ ai„iai>rn'Yimin)Yi'm'in + 0)). 

Im I'm' 

From isotropy we can define Ci from 

{airnai'm') = ^U'Smm'Cl (18) 

and therefore 

21 + 1 
47r 



uj{e) =Y,CiYira{n)Yi„,{h + e) ^Y,Ci^—Pi{cose), (19) 



Ir. 
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where the addition theorem was used in the last equahty. 

The Ci 's are the variance of the distribution for the coefficients aim i assumed gaussian: 

P{aim)daim = — ==e daim (20) 

The covariance matrix is given by: 

Cow{uj{e)uj{e')) = {uj{e)uj{e')) = (2f) 

^ Com{GiCv)^ -^Pii^os 9)Pi- (cos e') 

w 

For full sky the Cz's are statistically independent. For partial sky coverage one can use an averaged Ci over a band 
of width A/ = 10 and these will be independent to a good approximation. We will assume that 

Cov(QQO = Var(Q)<5„, (22) 

and assuming a gaussian distribution for the likelihood function of Ci [23| , correcting for partial sky coverage f^^y and 
including shot-noisc error one has: 

Var(CO=.^^^^(Q-M/.)^ (23) 

where n = N/AQ is the average number of galaxies per unit solid angle. 
Therefore the covariance matrix can be computed as: 

2 2/4-1 

Cov(w(0)c^(0')) = ^ E -7T^Pi{cos9)Pi{cos9') {Ci + 1/nf (24) 

/sky J l^TTj 

We can estimate the Q's from a model for the power spectrum by first performing a Fourier transform in the 3-d 
density field in ea. p^ 

/j3iL 
^5(fc,z)e^'=-"'-W (25) 



and use the identity 



to compare with eq. (|16p and find 



^ik-nr ^ 4^ ^ ^l^^ ^f^^Y^m (fc) (26) 



/d^k 
j^Sik,z)AWji{kr)Yi:^{k) (27) 

Defining 

^{k)^ j dz <i>,{z)D{z)3i{kr{z))b{z) (28) 
with (5(fc,z) = D{z)6{k,{)) allows us to write: 

CI = (la^n = ljdk k^P^Lik) (k) (29) 

At this point one must also introduce the effects of redshift distortions, as done for the angular correlation function. 
We will adopt the prescription given in [U, H^l and add to the function ^'/(A;) another term ^'[(fc) that incorporates 
the redshift distortion given by: 

*;■'■(&) = / dz p{z)cj,,{z)D{z) 



{2l + 3){2l - 1) 



jiikr) (30) 



l{l-l) . ([+1)0 + 2). 



{2l-l){2l + iy ' ' 2l + l)(2l + 3y 

In Fig. ([2]) we show the impact of the photo-z error on the angular power spectrum for = 0.03,0.05 and 0.10. 
One can see how the power decreases with increasing photo-z error, as expected. We also have checked the numerical 
agreement of the angular correlation function computed from cq. (jl9p with the sum on spherical harmonics up to 
I = 1000 (to ensure numerical accuracy) with the definition given in cq.([T]). 
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FIG. 2: Photo-z error impact on the angular power spectrum for the redshift shell 0.95 < z < 1.00. The solid line correspond 
to no errors in the redshift determination, dashed line with = 0.03, dot-dashed line ~ 0.05 and dotted line 5a = 0.10. 



B. Covariance matrix including different shells 



The various redshift bins of the survey are not independent due to the errors in the photometric redshift. Of course 
the main correlations will occur for adjacent bins. Hence we need covariance matrices for 2 different redshifts i and j 



9; I 1 2/' -I- 1 

V cov(QC/, ) ^ -—Pi (cos e,,)Pi, (cos e„ 

— ' 47r 47r 



(31) 



where in this case, since the shot noise between shells are uncorrelated, we write 

2 / n2 



Cov(Q'C/,) = 



(c;^^ +l/n.(5,j) ill. 



(2/ + 1)./.., 

Therefore, in the general multibin case the full covariance matrix can be computed as: 



J sky ^ 



2/ + 1 

— jP;(cos e„) 



(47r)2 

P,(cos0™)(c;-'+lM(5, 



with 



Ci 



(32) 



(33) 



(34) 



In Fig([3]) it is shown the cross-correlation angular power spectrum C/'"' for the redshift shell 0.95 < z < 1.00 with 
itself (auto-correlation) and its 3 nearest neighbours on one side. One can sec that the cross-correlation for the nearest 
2 shells is enough to capture the relevant effect. Hence in the following we will include the cross-correlation only with 
the neighbouring 4 shells, 2 on each side. 



ERROR FORECASTS FOR COSMOLOGICAL PARAMETERS 



We will use the full shape of the angular correlation function uj{9) in our analysis to derive forecasts for the errors 
in the determination of cosmological parameters using a Fisher matrix approach. 
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FIG. 3: Cross-correlation angular power spectrum for the redshift shell 0.95 < z < 1.00 with its 3 next neighbours. The solid 
line correspond to auto-correlation. The photo-z error was set to 5c, = 0.03. 



The Fisher matrix approach ag; 
distribution of the parameters [23 



)roximates the full likelihood function of the parameters by a multivariate Gaussian 
M: 



C{{p}) (X cxp 



1 



that is, the Fisher matrix Fajj is the inverse of the covariancc matrix for the parameters: 



92 log C{{p}) 



dpadpii 



(35) 



(36) 



{p} 



where {p\ denotes the best-fit parameters. We do not know a priori what the best-fit parameters are, since this 
involves a complicated search in the space with dimensions given by the total number of parameters. We will compare 
in the Appendix a few results from Fisher matrix with brute force computation of the likelihood function. 

The most useful application of the Fisher matrix approach is to study forecasts of future experiments, determining 
the precision with which the parameters can be measured with respect to fiducial values. These fiducial values replace 
the unknown best fit parameters. 

The Fisher matrix in this case is given by: 



^ap = ^ [C \ 



dpa 



dpp 



(37) 



dpa dpp 



where C^-'jn = Cov(a;*(0„)w-'(0„i)) is the covariance matrix derived in the previous section and the derivatives are 
calculated at the fiducial set of parameters {p}. In general the first term dominates over the second one (33j and our 
results will be obtained using only the first term. This is a conservative result and further below we will show that 
the second term indeed has a small impact on the forecast. 



VI. RESULTS 



We have divided a DES-like survey into 20 redshift shells with width Az = 0.05 in the range 0.4 < z < 1.4. We 
justify this value of Az below. The angular binning varies with the redshift shell in order to select angles around the 
BAO peak position. The redshift and angular binnings are shown in Table (|T| for the case of roughly 60 < r < 160 
/i~^Mpc. The full correlation matrix in our case turns out to be a 445 x 445 matrix. 
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bins 


Dili WlLlLli 




0.40 < z < 0.45 


3 - 8 


26 


0.2 


60 < r < 160 


0.45 < z < 0.50 


3 - 8 


26 


0.2 


67 < r < 180 


0.50 < z < 0.55 


2.5 - 6.1 


25 


0.15 


61 < r < 149 


0.55 < z < 0.60 


2.5 - 6.1 


25 


0.15 


66 < r < 162 


0.60 < z < 0.65 


2 - 5.6 


25 


0.15 


57 < r < 159 


0.65 < z < 0.70 


2 - 5.6 


25 


0.15 


61 < r < 170 


0.70 < 2 < 0.75 


1.7 - 5.0 


23 


0.15 


55 < r < 162 


0.75 < a < 0.80 


1.7 - 5.0 


23 


0.15 


58 < r < 171 


0.80 < z < 0.85 


1.5 - 4.5 


21 


0.15 


54 < r < 161 


0.85 < z < 0.90 


1.5 - 4.5 


21 


0.15 


56 < r < 169 


0.90 < 2: < 0.95 


1.5 - 4.5 


21 


0.15 


59 < r < 177 


0.95 < z < 1.00 


1.5 - 4.0 


21 


0.15 


61 < r < 183 


1.00 < 2 < 1.05 


1.3 - 3.7 


21 


0.12 


55 < r < 157 


1.05 < 2 < 1.10 


1.3 - 3.7 


21 


0.12 


57 < r < 163 


1.10 < 2 < 1.15 


1.3 - 3.7 


21 


0.12 


59 < r < 168 


1.15 < 2: < 1.20 


1.2 - 3.1 


20 


0.1 


56 < r < 145 


1.20 < 2 < 1.25 


1.2 - 3.1 


20 


0.1 


58 < r < 150 


1.25 < 2 < 1.30 


1.2 - 3.1 


20 


0.1 


60 < r < 154 


1.30 < 2 < 1.35 


1.2 - 3.1 


20 


0.1 


61< r < 158 


1.35 < 2 < 1.40 


1.2 - 3.1 


20 


0.1 


63 < r < 162 



TABLE I: Binning in redshift and angles (in degrees) for roug hly 60 < r < 160 h'^Mpc 



The Ci 's were computed in a suite of c codes based on the free software Gnu Scientific Library [26[ . The evaluation 
of the integrals of the rapidly oscillating spherical Bessel functions are particularly demanding. We compute the Ci 's 
up to / = 1000 in order to obtain our results. 

The reduced covariance matrix for one bin (1.00 < z < 1.05), defined by — C„m / -\/C„n Cmm is shown in 

Fig.(|4]). In this case we used f^^y = 0.125 and n = 15/arcmin^. It can be seen that the errors in different angular 




I" 

■ o.g 

- 0.8 

- 0.7 

- 0.6 
0.5 
0.4 
0.3 
0.2 

—1 0.1 

1.5 2.0 2.5 3.0 3.5 
9{degrees) 



FIG. 4: Reduced covariance matrix for the bin 1.00 < z < 1.05. 

bins are highly correlated resulting in a quasi-singular covariance matrix. The magnitude of the first and second 
off-diagonal elements is typically 0.9 and 0.8. This is one of the major disadvantages of using the ACF to derive 
cosmological parameters. 

In order to compute the Fisher matrix we invert the covariance matrix using singular value decomposition [27| . 
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A. Redshift shell width and shot noise 



The narrower the redshift shell is, the larger signal is obtained for the angular correlation function, as illustrated in 
Fig-®. However, one must be concerned with the impact of shot- noise on the forecast of cosmological parameters, 
which we analyze in this subsection. 



Az=0.01 




0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 

9{degree) 



FIG. 5: Angular correlation function for a redshift bin centers at z — 1.1 for 4 different bin widths Az. 

In Fig. © we plot the la error forecast on the determination of the equation of state w as a function of the shell 
width when marginalizing over the parameters ilcdm, and bias, with the other parameters fixed at their fiducial 
values, and with a photo-z error given by cr^ = 0.03(1 + z). These forecasts were obtained for redshift shells centered 
at 3 mean redshift with different shell widths Az. As expected, one can see that for very narrow redshift shells the 
shot-noise error becomes dominant even for the DES expected number of galaxies, which we assume to be Ngai = 300 
million. 

The best value for the shell width is around Az ^ 0.05 (a similar result was found by [2^, where they found that 
the best shell width for the the redshift space distortion parameter is ~ 0.07). The blue dotted line shows the error 
forecast on w when the shot-noise term is neglected. As expected if shot-noise is absent, the result improves as the 
width gets smaller because the angular correlation has more signal for thinner shells. 

Another interesting feature of Fig. ([6]) is the flatness of the forecast error on w for Az ^ 0.05. One naively would 
expect that the constraint should degrade for very wide shells due to the washing out of the spatial correlation 



function when projecting on the sphere. However, as pointed out by Simpson et al [21], the corresponding top-hat 
shell that simulates a true redshift distribution, denoted by Azth , can be obtained from the convolution of the top- hat 
shell in photometric redshift space Az with the probability distribution for the photo-z error eq. ()12[) . which can be 
approximated by 



Azth = ^ Az"^ + 12^2 . (38) 

From this relation we find for Az = 0.01, 0.05, 0.10 and 0.15 the following top-hat values for the shell at z = 1.1 
(green dashed line in figure[6]), respectively: 0.219, 0.224, 0.240 and 0.265. Therefore, only for very wide redshift shells 
the true (spectroscopic) top-hat shell will have an impact in the cosmological parameters due to projection effects 
compared to very thin shells. Hence in what follows we set Az = 0.05 through the entire redshift range. 



B. Forecasts from the ACF 



The sensitivity to the cosmological parameters for each shell increases with redshift, as exemplified in Fig.©, where 
we show the errors in the equation of state obtained from each shell independently. This can also be seen in Fig.([8|), 
where we show the forecasts on the equation of state w and the dark matter density fJcrfm obtained from 3 independent 
redshift shells, 0.40 < z < 0.45, 0.85 < z < 0.90 and 1.35 < z < 1.40. We have marginalized over erg and bias (we 
allow for a free bias parameter for each shell) but the other parameters are held at their fiducial values. One notices 
that the forecasts are better for shells at larger redshifts and that the ellipses rotate slightly from one shell to the next, 
a behaviour which helps to break some of the degeneracies. The forecast errors are = 0.071 and oa^^,^ = 0.036. 
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FIG. 6; The lines show the expected error on lo as a function of the width of the shells centered at 3 mean redshifts. For all 
of them we are assuming Ngai = 300 millions of galaxies and = 0.03(1 + z). The red dot-dashed line is centered a.t z — 0.9, 
the dashed green for z = 1.1 and the solid blue line for z = 1.3. The dotted blue line is the same as the last shell but now 
neglecting the contribution from the shot-noise. 
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FIG. 7: Forecasts for l-cr errors on the equation of state w as a, function of redshift with marginalized bias, as and Qcdm- 
Optimistic case with Sa = 0.03 and 20 < r < 200 /i^^Mpc (sohd green line), fiducial case with (5o- ~ 0.03 and 60 < r < 160 
ft~^Mpc (dashed blue line) and conservative case (5a = 0.05 and 60 < r < 160 ft^^Mpc (dot-dashed red line). 

We also checked the robustness of this result if the other parameters are marginalized instead of fixed to their 
fiducial values. This is shown in Fig.®, where the previous result with all parameters other than J7cdm and w fixed 
at their fiducial values is shown in the solid ellipse. As expected, marginalizing over all other parameters without any 
priors (dot-dashed ellipse) degrades the forecasts significantly. The errors on the parameters without any priors are: 
cTu, = 0.34, (Tned„ = 0.047, CTOf, = 0.022, c/i = 0.26, cr„^ = 0.31 and a^^ = 0.77. However, introducing priors only on h 
(from HST) and il.b from WMAP7 (dashed ellipse) again reduces the errors of the forecasts. 

As explained above, our results will be conservative since we are neglecting an extra term in the Fisher matrix. To 
assess the small impact of the neglected term, we plot in Fig. ([T0| the analysis for all redshift bins without correlations 
among bins with and without the contribution from the neglected term. We sec that, as expected, including the 
second term leads to a slightly more restrictive forecast. 

We now proceed to our complete analysis. We consider 26 parameters, namely n^, as, h, 57f,, ilcdm,w, ^i, ^2, • • • , 
where bi are the bias for the i*'' redshift shell. We include correlations for the adjacent 2 redshift bins, i.e. for each 
bin we take correlations with 4 bins. In order to show the impact of including these correlations we show in Fig. pTj) 
the 1(7 constraint on ficdm and w with all other parameters marginalized without priors. As expected, it can be seen 
that the error is underestimated when the correlations are not taken into account. 

In Fig. (|12l) we plot the la ellipses for the parameters Us, fib, flcdm and w, marginalizing over (Tg and bias. When 
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FIG. 8: Forecasts for Qcdm and w from the angular correlation function of three redshift shells, 0.40 < z < 0.45 (blue dot-dashed 
line), 0.85 < z < 0.90 (green dashed line) and 1.35 < z < 1.40 (red solid line). The yellow ellipse is the result of combining the 
3 shells. 




FIG. 9: Forecasts for Slcdm and w from the angular correlation function of five redshift shells. The red solid ellipse is the result 
with all parameters fixed at their fiducial values; blue dot-dashed ellipse is the result with all parameters marginalized without 
priors; green dashed ellipse is the result with priors on h and fi;,. 



two of these parameters are plotted all the others are marginalized with no priors. We study three scenarios: optimistie 
(20 < r < 200 /i-iMpc and S„ = 0.03), fiducial (60 < r < 160 h'^Mpc and = 0.03) and pessimistic (60 < r < 160 
h~^Mpc and So- = 0.05). One can see that a larger photo-z error degrades substantially the possible forecasts. This is 
expected since for larger photo-z errors the signal coming from the spatial correlation function becomes weaker. The 
forecast error for w goes from ^ 0.3 with = 0.03 to ^ 0.5 for da- ~ 0.05. 

One can also notice from Fig. (|12p that the constraint in w is not very sensitive to the scales probed. Nevertheless, 
the other parameters ficrfm, ^^6, and h are significantly dependent on the scales probed, especially and Ug. This 
shows that those cosmological parameters change the ACF at all scales and not only around the BAO feature like w. 

Finally we investigate the role of priors in the forecasts in Fig. p^ . In particular, wc use the HST prior on h [2^ 
and WMAP7 priors from the acoustic scale Ia , the redshift of decoupling and the shift parameter R (table 10 of 
Komatsu et al [3) and also imposed a prior in rig from WMAP7. The WMAP7 covariance matrix in I a, and R 
was then converted to our set of cosmological parameters using a method described in Mukherjee et al [s^l and Wang 
et al 

We find the l-cr marginalized error on a given cosmological parameter p via 

= Vt^^'lpp- (39) 
In our optimistic best case scenario one would be able to obtain — 0.2 and cro^^^, = 0.03 without priors and 
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FIG. 10: Forecasts for Qcdm and w from the angular correlation function for all redshift shells without correlation among shells. 
Contribution from the neglected term (dashed line), contribution from the dominant term (dotted line) and result from both 
terms (solid line). 



0.28 



0.26 

0.24 

g 

0.22 
0.20 
0.18 



-1.4 -1.3 -1.2 -1.1 -1.0 -0.9 -0.8 -0.7 -0.6 
W 



FIG. 11: Forecasts for ficdm and w from the angular correlation function with all redshift shells included. All other parameters 
are marginalized with no priors and Scr ~ 0.03. Dashed line is the result with diagonal correlation matrix (only auto-correlations) 
and full line is the result taking into account next 2 neighbour bins. 



a-ui = 0.08 and <Jn^am ~ 0.009 with priors. 



C. Impact of fiducial bias 

The analysis presented until now was performed assuming a fiducial galaxy bias 6 = 2 for each redshift shell. This 
is the usual value for catalogues of luminous red galaxies (LRG). However, it is more likely that the galaxy bias will 
have lower values for the complete galaxy set. In this subsection we show how our results change if the galaxies in 
the catalogue are assumed to be unbiased on average with 6 = 1, as done in [28| . 

As shown in [s^ , the statistical reach of a survey can be assessed from its effective volume, which is proportional to 
b^. Hence, it is expected that the constrain on the cosmological parameters will degrade if the fiducial bias decreases. 
In order to evaluate the impact of the bias on the cosmological parameters, we have performed the Fisher matrix 
analysis with a fiducial bias b ^ 1 and compared to the results arising from 5=2. 

In Fig.p^. we show the forecast for our optimistic scenario. In this case the constrain on w degrades by 17%, with 
respect to 6 = 2, and the expected constrain on w is cr^ = 0.24. This analysis shows that the fiducial galaxy bias does 
have a minor impact in the dark energy constraints. 
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FIG. 12: Icr forecasts for ns,h,Qb,Q,cdm and w, marginalizing over as and bias, for three scenarios: optimistic (solid line) 
(20 < r < 200 Mpc and 5^ = 0.03), fiducial (dashed hne) (60 < r < 160 Mpc and 5^ = 0.03) and pessimistic (dooted line) 
(60 < r < 160 Mpc and 5^ = 0.05). 




FIG. 13: Qcdm X TO 1(7 forecasts in the optimistic scenario in the case of no priors (dashed line), and adding priors from WMAP7 
and HST as explained in the text (solid line). 



VII. SUMMARY AND CONCLUSIONS 



In this paper we have investigated the forecasts on cosmological parameters that can be obtained from a study of the 
full shape of the 2-point angular correlation function that will be measured by future large-scale photometric surveys, 
such as the Dark Energy Survey project. The angular correlation function was modelled taking into account redshift 
space distortion, nonlinear corrections to the power spectrum, bias and gaussian photo-z errors. After a proper binning 
in angle and redshift, the Fisher information matrix was constructed from a covariance matrix that includes correlation 
of a given redshift bin with the 4 nearest bins. We considered 26 parameters: ns,as^ h, fib, i^cdrm 'w,bi^b2, ■ ■ ■ , &20 with 
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FIG. 14: Icr forecats for w and ilcdm for the optimistic scenario with two fiducial bias: b = 2 (blue solid line) and 6=1 (red 
dashed line). 

different bias for each shell. 

We used the 26 x 26 Fisher matrix to obtain Icr ellipses for the cosmological parameters around the adopted fiducial 
values obtained from the central values of WMAP7 measurements. In all our results we marginalize over the bias 
parameters and erg. After showing results for several combinations of parameters in Fig. (jl2p . we concentrate on 
forecasts on the flcdm x w plane. We have shown the improvements of using priors from other experiments such as 
WMAP7 and HST in comparison to a simple marginalization of parameters. Finally, we showed that the effect of 
neg lecting the second term in the Fisher matrix cq. (|38p is in fact unimportant, as usually assumed in the literature 

In the appendix we make a comparison between the Fisher matrix approximated methodology and a full likelihood 
search in parameter space which validates the Fisher matrix approach. 

We find that under these assumptions the Dark Energy Survey should be able to constrain the dark energy equation 
of state parameter w and the cold dark matter density ^cdm with a precison of 20% and 13% respectively from the 
full shape of the angular correlation function alone. When combined with priors from other observations the precision 
in the determination of these parameters increase to 8% and 4% respectively. 

For comparison, the WiggleZ Dark Energy Survey has recently obtained w = —1.6^0 7 from BAO data alone, which 
is improved to w — — 0.982lQ'i[g4 when WMPA7 distance priors are included [l^l, which amounts to a precision of 
approximately 40% (15% with priors). Therefore, measurements of the angular correlation function from a DES-like 
survey will improve the constrain on w by a factor of C(2), in the case where only galaxy clustering measurements 
are used. It is also expected that the constraints we have found will improve as the redshift range, area and number 
of galaxies improve with future photometric surveys. 

In conclusion, our results indicate that an analysis of the full shape of the 2-point angular correlation function can in 
fact bring an important contribution to the determination of cosmological parameters in future photometric redshift 
surveys. This contribution is of course complementary to the other observables that will be used in a DES-like survey, 
including cluster number counts, shear from weak lensing and distances from SNIa. One challange that remains is to 
combine our results with these other probes in a consistent manner, including the proper correlations. 

Appendix 

In this Appendix we check the reliability of the Fisher matrix forecast by computing the likelihood function on a 
grid of values for the parameters. Here we concentrate on a single redshift shell, 0.70 < z < 0.75, and examine 4 free 
parameters, namely w, J^cdm, b and cts, with a grid constructed around their fiducial values, see Table Once 
the grid is created it is fairly easy to search the likelihood function and compare it with the Fisher matrix approach. 

In Fig. (|15p we compare the likelihood for w when marginalizing over the other 3 parameters. In this case, the 
Fisher matrix prediction is in good agreement with the results from the grid. The Fisher matrix forecasts an error of 
o'wifisher) = 0.23, exactly the same value obtained from the grid. It is interesting to note that the true distribution 
has a small skewness towards higher values, showing that it is not a perfect gaussian distribution, but in overall both 
distributions are in fair agreement. 
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Parameter 


lower bound 


upper bound 


grid 


w 


-2.0 


0.0 


80 


^CDM 


0.05 


0.55 


100 




0.1 


2.1 


100 


b 


0.05 


5.55 


100 



TABLE II: Grid parameters. 



For flcDM we find the same trend as for w, but here the difference between Fisher matrix and the grid likelihood is 
slightly larger, as shown in Fig. (jl6p . In this case the Vic dm has a more pronounced non-gaussianity, with a skewness 
towards higher values. However, the prediction for the standard deviation from Fisher matrix is again in fairly good 
agreement with the grid, with ucDMifisher) = 0.088 and crcDM^gTid) = 0.064. In this case the Fisher matrix 
method overestimates the error on Vic dm due to the non-gaussianity of the likelihood. 

This simple comparison shows that the Fisher matrix gives reliable predictions for the forecasts of the cosmological 
parameters w and ^cdm that we obtained in this work. 



2.0 




FIG. 15: The actual distribution for w estimated with the grid when the other 3 free parameters are marginalized (solid line) 
compared with the Gaussian distribution predicted by the Fisher matrix (dashed line). 




FIG. 16: The actual distribution for ^cdm estimated with the grid when the other 3 free parameters are marginalized (solid 
line) compared with the Gaussian distribution predicted by the Fisher matrix (dashed line). 
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